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ABSTRACT 



The steam generator of a DLG-9 Class Destroyer is studied and 
simulated by means of the Digital Simulation Language DSL/360. 
Parameter plane studies are made for each one of the principal subloops 
of the control system and determination of optimal controller settings 
is attempted . 

The entire steam generator and control system are simulated and 
the responses to small disturbances are analyzed. 

Finally a linear interpolation is attempted for the nonlinear 
transfer functions, and the responses compared with data from an 



actual DLG-9 Test Boiler. 
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I. INTRODUCTION 



A. DISCUSSION 

The Naval Boiler and Turbine Laboratory has conducted extensive 
dynamic and frequency response tests on a Babcock and Wilcox 
Company test boiler of the DLG-9 class, and through the use of the 
"pulse test" technique described in Ref. 1 a mathematical model was 
obtained of the steam generator system with its controls and auxiliary 
equipment . 

This mathematical model is reported in NBTL PROJECT B-502-III 
and Figure 1, taken from that report is the block diagram of the model 
with its associated transfer functions. It is the objective of this 
thesis to simulate the steam generator using the Digital Simulation 
Language DSL/3 60 [2], with the transfer functions as reported in 
Ref. 1 and to establish the validity of the model; also using parameter 
plane techniques [3] to perform the analysis of each of the three 
principal sub-loops (air flow, oil flow and water level control systems) 
in the decoupled state. An attempt to improve the response of these 
loops will be made by changing the controllers gains; however, since 
the specifications for the performance of the loops in the decoupled 
state are not known, those specifications will be defined with the 
criterion of making the step response as fast as possible with little 
or no overshoot. A simulation of the complete steam generator will be 
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made with the parameter settings as chosen from above specifications 



and the response compared with those for the parameter settings as 
given in Ref. 1 . 

Finally a linear interpolation will be attempted for the nonlinear 
devices and responses of the simulation using the linear interpolation 
will be compared with actual responses of the boiler as reported in 
Ref. 4. 

B. DESCRIPTION 

The DLG-9 boiler utilizes the Bailey Combustion Control System, 
of which a general description will be given in this section, and a more 
detailed description will be given in later sections. 

The primary functions of the automatic boiler control system are: 
to control the firing rate of the boilers to maintain a constant super- 
heater outlet pressure under all conditions of boiler steam load; to 
proportion the combustion air (quantity of air needed to burn the fuel oil 
furnished to the boiler) and fuel in the proper ratio to maintain the 
optimum combustion efficiency; to reduce the firing rate to a stable 
minimum in the event of instrument failure; finally, through the feed- 
water control, to regulate the flow of feedwater to the boiler drum to 
maintain a normal boiler drum level. 

The amount of fuel burned, the quantity of air required to burn 
that fuel and the amount of feedwater needed to maintain the boiler 
drum level, are determined by the amount of steam being used and the 
desired pressure to be maintained. 
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The amount of steam output demanded of the boiler, Gs, is 
measured, through its entire range (0 to 210000 Ibs/hr) , by the steam 
flow transmitter which develops a pneumatic signal, Pgs, (3 to 2 7 psi) 
proportional to the steam flow, which is sent, as demand signal, to 
the steam pressure controller and the water level control loop. 

The desired superheater outlet pressure to be maintained is 
1200 psi. This pressure is measured by the superheater steam pressure 
transmitter which develops a pneumatic signal, Pp, from 3 to 27 psi, 
in direct proportion to a range of superheater pressure variation from 
900 psi to 1500 psi; therefore to 1200 psi of superheater pressure 
corresponds a signal of 15 psi which is compared to a reference or set 
point, Prp, of 15 psi also, so as to have a signal. Pep, that is a 
measure of the deviation of superheater pressure from 1200 psi. This 
signal is also sent to the steam pressure controller whose output depends 
on the two signals mentioned. Pep and Pgs. The output of the con- 
troller is sent to the fuel oil control loop to modify the firing rate 
according to the demand, and also is sent to air control loop to modify 
the amount of air flow in accordance to the change in the firing rate 
which is going to take place. 
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Steam Flow 




FIGURE 1 

SCHEMATIC DIAGRAM OF BOILER CONTROLS 
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TABLE I 



SET 



SPCl 



SPC2 



AFC 



RABR 



FDBA 



FDB 



AFT 



CR 



TRANSFER FUNCTIONS CORRESPONDING TO FIGURE 1 

Steam Flow Transmitter: 

(1 . 14xl0"^)/(l + 1.36S + 0.052S^) 

Steam Pressure Controller: 

(2.5) /(l + 0.034S) 

Steam Pressure Controller: 

(1 .0)/(l + 0.038S + 0.004S^) 

Air Flow Controller: 

(0.944 + 0.106/S)/(1 + 0.04S + 0.003S^) 
Rate Action Booster Relay: 

(1 + 0.895)7(1 + 0.079S + O.OIS^) 

Forced Draft Blower Actuators #1 and #2: 

Crsng: (512/(1+S))/(1 + 0.28S + 0.5S^) 

90% FP: (950/(l+S))/(l + 0.28S + 0.5S^) 

Forced Draft Blowers #1 and #2: 

Crsng: (6.83x10 ^)/(l + 10. OS) 

90% FP: (5.50xl0”^)/(l + 6.35S) 

Air Flow Transmitter: 

(0.254exp(-0.45S))/(l + 0.4S) 
Calibrating Relay: 

(1.05) 7(1 + 0.034S) 
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EAB 



Excess Air Bias: 



SSR 



OFC 



FOCV 



ROFT 



SOFT 



OFR 



SHPT 



SH 



WLC 



FWRV 



Amount of excess air to be added or subtracted 
to produce smoke or to blow tubes. 

Signal Selector Relay: 

Transmits the smaller value of Pbq and Pm 
Oil Flow Controller: 

(2.82 + 0.392/S)/(l + 0.32S) 

Fuel Oil Flow Control Valve: 

795.0 

Return Oil Flow Transmitter: 

(1.2xl0"^)/(l + 0.12S + 0.33S^) 

Supply Oil Flow Transmitter: 

Same as ROFT 
Oil Flow Relay: 

(1 .05)/(l + 0.034S) 

Superheater Steam Pressure Transmitter: 
(4.0xl0”^)/(l + 0.012S + 0.02S^) 

Superheater: 

-4 

Crsng: 3.0x10 
90% FP: 8.2xl0”'* 

Water Level Controller: 

(1.26 + 0. 1/S)/(1 + 0.54S + 0.8S^) 

Feedwater Regulating Valve: 

(3.12exp(-0.25S))/(l + 0.68S + 1.26S^) 
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WFT 



Water Flow Transmitter: 



F 



DR 



WLT 



SP/SF 



SP/OF 



WL/S 



WL/S 



SP/WF 



(1.14x10 ^)/(l + 1.36S + 0.052S^) 

Filter: 

(1.0)/(1 + 4.5S) 

Steam Flow-Water Flow Differential Relay: 
(0.50)/(0.034S) 

Water Level Transmitter: 

(1.0)/(1 + 0.6S + 0.25S^) 

Steam Pressure / Steam Flow: 

Crsng: (1 . eSKlo”"^) (1+35 . 70S)/(S(l+22 . 5S) (1+0 . 71S)) 

90% FP: (2 .5xl0”^)(l+0.504S)/(S(l+2.23S++ 2.53S^)) 
/ 

Steam Pressure / Oil Flow: 

Crsng: (2 . 83xlo”^) (1+35 . 6S)/(S(l+22 . 5S) (1+0 . 57S)) 
90% FP: (3 . 18xl0”^) (1+4 . 50S)/(S(l+2 . 50S) (1+0 . 71S)) 
Water Level / Steam Flow: 

Crsng: (1.12x10 ^)/S 
90% FP: (1.26xl0"^)/S 
Water Level / Steam Flow: 

Crsng: (5 . 62xl0"^)/((l+3 . 18S) (1+0 .48S+0 . 63S^)) 

90% FP: (4.04xl0~^)/((l+1.80S)(l+0.21S+0. 13S^)) 
Steam Pressure / Water Flow: 

Crsng: (0 .389xl0"^)/(S(l +22.5S)) 

90% FP: (0. 100xl0"'^)/(S(l + 6.0S + 100. OS^)) 
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WL/W 



WL/0 



Water Level / Water Flow: 

Crsng: (1 . 12xlO"®)/(S(l + 0.852S + 2.02S^)) 

90% FP: (1 .26xl0”®)/(S{l + 0.320S + 0.63S^)) 

Water Level / Oil Flow; 

Crsng: (2 . 09x1 0'^^expC-0 . 8S)/((l+0 .453) (1+1 . 9S+2 . 53^)) 
90% FP; (1 .9xl0~^exp(-0.363)/((l+2 .03) (1+2 .03+4.03^)) 



20 



II, AIR FLOW CONTROL SYSTEM 



A. DESCRIPTION 

The Air Flow Control System is shown in block diagram in 
Figure 2; it is a typical feedback system with two principal non- 
linearities, the steam control valves and the forced draft blowers. 

The system consists of the following elements: 

1 . The Controller 

The Controller is a Bailey Meter Company "Mini-line" 
Standatrol, proportional plus reset controller; it has as input the 
loading pressure, Pe, that represents the difference between the 
master demand signal. Pm, from the steam pressure controller, and 
the signal from the air flow measurement. The controller is calibrated 
to maintain a constant output pressure as long as Pe is zero or when 
the air flow measurement is equal to the demand signal. 

2 . Rate Action Booster Relay 

A Bailey Meter Company Rate Booster Relay, which is a pro- 
portional plus rate action controller, is introduced to provide 
increased phase margin in the blower control loop, thus allowing 
maximum loop gain. It is driven by the air flow controller, and if its 
input changes, the rate action booster relay will commence to in- 
crease or decrease its output at the rate set by the proportional band 
and reset settings until the air flow again balances the demand signal. 
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FIGURE 2 

BLOCK DIAGRAM. AIR FLOW CONTROL SYSTEM 




The output of this controller serves to position the forced draft blower 



throttles and the air flow control dampers. 

3 . Forced Draft Blower Actuator 

Linkage and throttle valve, Bailey Meter Company Model 
AC-44 Control Drive and Mason-Neslan Control Valve, constitutes the 
steam regulating system for the main forced draft blowers. The 
throttle valve which admits the steam, Gb, to the forced draft blower’s 
turbine is positioned according to the pneumatic signal, Pq, from the 
booster relay. The dynamics of the actuator are assumed to be linear 
for all blower speeds but the gain varies as the blower's speed is 
varied . This nonlinearity is due to the characteristics of the V-ported 
blower steam control valve. 

4 . Main Forced Draft Blower 

Results obtained (NBTL Report) by generating a triangular 
pulse in the pneumatic signal to the actuators, indicates that the 
blowers can be represented by a first-order lag characterized by a 
natural frequency increasing with load, and a gain decreasing with 
load. This nonlinearity is due to the driving torque and retarding 
torque being nonlinear functions of turbine steam flow and blower 
speed . 

5 . Air Flow Transmitter 

This is a Bailey Meter Company Type CJ-20 Differential 
Transmitter with Type KC-16 square root converter - transmitter. 
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This instrument in the feedback circuit of the air flow control loop, 
measures the pressure differential of the air flow across an orifice, 
extracts the square root of the differential pressure, and develops a 
pneumatic loading pressure that is proportional to the flow of air. 

6 . Air Flow Calibrating Relay 

This relay, operating on the air flow transmitter output 
signal, is provided to permit scaling of the air flow open-loop gain in 
order to obtain optimum air/fuel ratio. 

B . ANALYSIS AND SIMULATION 

The Air Flow system was simulated using DSL, and using the 
transfer functions as given in Figure 2. The transient response to a 
unit step input is shown in Figure 3 for 90% Full Power condition and 
in Figure 4 for Cruising condition. As noted from Figure 3, there is 
some oscillation in the transient, and the settling time is about 30 
seconds; the Cruising condition response is slower, about 40 seconds 
settling time. 

It is desired to "optimize" the response of the loop, by adjusting 
the coefficients of the Air Flow Controller; therefore the Controller is 
defined as: 

Alpha + Beta / s 

where Alpha and Beta are the variable coefficients, which in the NBTL 
report are given as: 

Alpha = 0.944 
Beta = 0.106 
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PERCENTAGE AIR FLOW 




FIGURE 3 . UNIT STEP RESPONSE 
AIR LOOP (907oFP). ALPHA = 0.944, BETA = 0.106 
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AIR FLOW (PERCENTAGE) 




FIGURE 4 . UNIT STEP RESPONSE 
AIR LOOP (CRUISING) . ALPHA = 0.944, BETA =0.106 
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By optimizing the loop response it is assumed that the desired response 
is a faster one, and in the case of the 90% FP a smoother response is 
desired . 

In order to calculate the characteristic equation for the loop a 
Fade approximation of the fourth order was made to represent the time 
delay shown in Figure 2; the Fade approximation in question is; 

, . 1680 - 840x + 180x^ - 20x^ + 4x^ 

exp(-x) = ^ — 

1680 + 840x + 180x + 20x + x 

where x = 0.045s. 

The resulting characteristic equation was of 15th order. Since 
there are two variable coefficients the best method of analysis is the 
parameter plane which is a graphical method that gives the following 
curves as function of the two parameters: 

Constant zeta curves (maps of radial lines in the s-plane) for 
the range of frequency Wn, specified by the user. 

Constant Wn curves (maps of circles centered at the origin of 
the s-plane) for a pre-programmed set of zeta values. 

Constant sigma curves, each one of which is the map of a 
specific point in the real axis of the s-plane. 

Constant zeta-Wn product curves. 

Figure 5 is the parameter plane plot for the 90% FP condition. 

In this graph, for clarity only, the constant zeta curves were drawn 
to show that some regions of the graph are covered by two sets of 
zeta lines; this is so because within the range of Wn of interest there 
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are two pairs of complex roots . This is more clearly seen in the Root 



Locus plot of Figure 6, made holding Beta constant at the given value 
of 0.106, and varying Alpha for 3 decades. The roots marked with a 
circle are at Alpha = 0.944 and those marked with a triangle are at 
Alpha = 0.60. It is seen that as Alpha decreases from 0.944 to 0.6 
the real roots closer to the origin become complex. This is seen in 
the parameter plane Figure 5 by moving the operating point ”B" 

(Alpha = 0 . 944 , Beta = 0 . 106) in a horizontal line to the left, i . e . , 
decreasing Alpha with Beta constant. At a point where Alpha = 0.66 
begin the second set of zeta lines, zeta = 0.99, which indicates that 
at this point the real roots closer to the origin are leaving the real 
axis . 

The outlined region in Figure 5 is expanded in Figure 7, which 
shows the constant zeta, sigma and Wn lines. 

Figure 8 is the parameter plane plot for the cruising condition. 
Point "B" in Figures 7 and 8 indicates the operating point with 
parameter settings as given by the NBTL report. As stated in the 
Introduction, it is desired to reduce the rise time and have little or 
no overshoot; in Figure 7 it is seen that the dominant roots are real 
and at sigma = -0.1 and sigma = -0.35, which determine the general 
shape of the response. The high frequency oscillation noted in 
Figure 3 is due to the effect of the complex roots at zeta = 0.158 and 
Wn = 1.31; to speed up the rise time, zeta should be decreased, but 
for the 90% FP condition the rise time is determined by the real roots 
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BETA 




FIGURE 7. PARAMETER PLANE 
AIR LOOP (90%FP) 
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BETA 




FIGURE 8. PARAMETER PLANE 
AIR LOOP (CRUISING) 
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that are dominant, therefore, by choosing an operating point like "C" 
(Alpha = 0.68, Beta = 0.12) where the dominant roots have become 
complex with zeta = 0.9 and Wn = 0.18, the zeta have been effectively 
lowered hence decreasing the rise time, ^Iso it is noted that the zeta 
for the second pair of complex roots have been increased from 
zeta = 0.157 to zeta = 0.178, a very small increment which is believed 
to diminish the oscillation noted in Figure 3. The step response for 
operating point "C" and for 90% FP is shown in Figure 9 where it is 
compared with the response for operating "B" to show that the rise 
time has been decreased. Examination of Figure 8 shows that point 
"C" as chosen above will also lower the zeta for this, the cruising 
condition, but in this case a small overshoot appears as shown in 
Figure 10. 
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AIR FLOW (PERCENTAGE) 




TIME (SECONDS) 

FIGURE 9. UNIT STEP RESPONSE 
AIR LOOP (90%FP). ALPHA = 0.68, BETA = 0.12 
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AIR FLOW (PERCENTAGE) 
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III. FUEL OIL CONTROL SYSTEM 



A. DESCRIPTION 

The fuel oil flow control system is shown in block diagram in 
Figure 11, it has as input the pneumatic signal, Pdf, as the demand 
index and the output is the fuel flow of the boiler, Gf. There is a 
constant input of 15900 Ibs/hr. of fuel oil supplied to the system at a 
constant pressure of 1000 psi; effectively the amount supplied to the 
atomizers is not exactly 15900 Ibs/hr. minus Gf, but less, due to 
back-pressure, since some pressure is needed to atomize the fuel oil; 
this is represented in the block diagram by making Gp equal to 15900 
minus 0.0445 Gf; Gp being the effective supply of fuel oil. The 
system is of the return flow burner type, which means that the control 
is in the return line, done by the fuel oil control valve; the oil flow 
control system is composed of the following elements; 

1 . The Controller 

This controller, as the one in the air flow control system, 
is a typical proportional plus reset "Mini-Line" Standatrol controller. 
The inputs to the controller are the pneumatic signals Pf, from the oil 
flow relay and Pdf from the fuel limiting selector relay; these two 
signals are balanced by the controller and if they differ an output 
signal, Pdr, is generated that is the control signal to the return fuel 
oil flow valve . 
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FIGURE 11 

BLOCK DIAGRAM. FUEL OIL FLOW CONTROL SYSTEM 





2 . Return Fuel Flow Control Valve 

This is a diaphragm-operated, V-ported valve, located in 
the fuel oil return line and controls the amount of fuel oil returned 
according to the input control signal. Note that the transfer function 
for this valve given in Figure 11 is different from the one shown in 
Figure 1; this discrepancy will be explained in the Analysis and Simu- 
lation part of this system. 

3 . Supply and Return Fuel Flow Transmitter 

These two transmitters are area meters that measure the 
supply fuel flow Gp (Ibs/hr) , and the return fuel flow Gr (Ibs/hr) , by 
means of a metering valve, and each one transmits a pneumatic signal, 
representative of the amount of fuel flow. 

4 . Oil Flow Relay 

This relay receives as input a signal which is obtained by 
subtracting the signal Pr, the output of the return fuel flow transmitter, 
from the signal Pp produced by the supply fuel flow transmitter. The 
relay develops an output, Pf, that is linearly proportional to the flow 
of fuel oil burned. 

B. ANALYSIS AND SIMULATION 

The oil flow control system when simulated on the digital 
computer using the transfer functions furnished by NBTL was unstable. 
The following block diagram reductions were made in order to find 
the characteristic equation of the loop and investigate the cause of 



the instability; 
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Gf 




where transfer functions for the different components of Figure 11 are 
defined by the letters inside the corresponding blocks. 

Moving the summing point 'a' to the input side of the trans- 
mitters, and reducing the two transmitters to one block since Grt is 
equal to Gst, the above block diagram is reduced to the following one: 
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It is noted that the input to the Gst block is Gf; therefore, the block 
diagram can be reduced to: 




In order to eliminate the implicit loop around the constant 0.0445, the 



following substitutions were made; 

Gf = Gp - Gr (1) 

Gp = 15900 - 0.0445Gf (2) 

Substituting (2) in (1) 

Gf = 15900 - 0.0445Gf - Gr (3) 

Gf = (15900 - Gr)/1. 0445 (4) 

And the block diagram was finally reduced to: 
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From this block diagram the characteristic equation was calculated as: 
GcGvGdGst - 1.0445 = 0 (5) 

When the numbers given by the report of NBTL for the indicated 
transfer functions were substituted in (5) , the negative sign persisted 
giving a characteristic equation of mixed signs which indicates roots 
in the right half plane; therefore, a change in sign was indicated in 
the fuel oil flow loop in order to have a potentially stable system. 

This change of sign was made in the return valve, since an increase 
in the demand signal produces an increase in the control signal, Pdr, 
and the return valve should produce a decrease in the amount of fuel 
oil returned in order to meet the demand of more oil to be burned. For 
this action to take place the sign of the return valve has to be negative. 

Making the indicated change in sign a parameter plane was 
plotted. Figure 12, for the resulting characteristic equation, and as 
before assuming the controller gains as variables, in this case: 

Gamma + Delta / s 

The values given for these parameters in the NBTL report are: 

Gamma = 2.82 
Delta = 0.392 

It is clearly seen in Figure 12 that for the values given above, 
operating point marked "A" in the figure, the system is unstable 
because the operating point is outside of the region covered by the 
zeta curves, meaning that the complex roots are in the right half part 
of the s-plane. To make the system stable, the gain Gamma could be 
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DELTA 




FIGURE 12. PARAMETER PLANE 
FUEL OIL LOOP 
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reduced to a value less than 2.0, keeping Delta constant at the given 
value. A simulation was done with Gamma set to 1.8, which showed that 
the system was stable for this parameter setting, as expected. The 
author of this thesis, thinking that such a gross instability as the one 
exhibited by this loop could not have been unnoticed by the authors of 
the NBTL report, through personal letter brought his findings to their 
attention. As a result, the transfer function for the control valve was 
changed as explained in Appendix A. 

Using, for the fuel oil return control valve, the transfer function 
indicated in Appendix A, a parameter plane was plotted. Figure 14. 

Again in this figure, the point marked "B" is for the parameters as 
given in NBTL report; the roots for this operating point are at zeta = 
0.145, Wn = 4.93 and sigma = 0.075. 

The fuel oil control system has a fast response, as seen on 
Figure 13, with some oscillation in the transient. If, as before, it 
is desired to eliminate this oscillation, the zeta of the complex roots 
should be increased, which is obtained by decreasing Gamma. But 
as Gamma is decreased, with Delta constant, Wn decreases and 
sigma does not change appreciably, and since the real root is so close 
to the origin, no significant change is obtained in the response. 
Therefore, it is necessary to move the real root farther from the origin, 
so as to make it less dominant, and let the variation in zeta have some 

noticeable effect. To increase sigma, Delta should be increased, so 



43 




FIGURE 13. UNIT STEP RESPONSE 
FUEL OIL LOOP. GAMMA = 2.82, DELTA = 0.392 
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GAMMA 

FIGURE 14 . PARAMETER PLANE 
FUEL OIL LOOP (CORRECTED DYNAMICS) 
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by choosing an operating point like "C" (Gamma - 2.0, Delta = 0.8) 
where zeta = 0.195, Wn = 4.94, and sigma =0.185, a unit step 
response is obtained, shown in Figure 15, that is somewhat faster 
than the original and less oscillatory, as was desired. 
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FIGURE 15. UNIT STEP RESPONSE 
FUEL OIL LOOP . GAMMA = 2.0, DELTA =0.80 
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IV. WATER LEVEL CONTROL SYSTEM 



A. DESCRIPTION 

The water level control system is a conventional three-element 
feedwater control system. It is designed so that the proportional plus 
reset control mode regulates feedwater flow rate while simultaneously 
maintaining boiler drum water level at the desired value. 

Figure 16 is the block diagram of the system. It has as demand 
index the signal Pgs that corresponds to the steam flow; it also has as 
input the signal PI that represents the water level of the boiler drum, 
and the reference signal Prl which is set to 15 psi as before. 

Note that the transfer functions of the controller, the feedwater 
regulating valve and the water flow transmitter in Figure 16 are dif- 
ferent from those given in Figure 1. The changes for the dynamics of 
the controller and water flow transmitter were made in accordance to 
Appendix B, where the thesis referenced in paragraph c of the appendix 
is Ref. 5 of this thesis. The correction for the gain of the feedwater 
regulating valve was made in accordance to Appendix B. 

The system consists of the following elements: 

1 . The Controller 

This controller, also a Standatrol, is a three-element 
controller which compares the signal Psw' from the steam flow-water 
flow differential relay, and the drum water level signal, PI, against 
the reference signal Prl and develops a pneumatic control signal, Pdw, 
sent to the feedwater regulating valve. 
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FIGURE 16 

BLOCK DIAGRAM. WATER LEVEL CONTROL SYSTEM 



Feedwater Regulating Valve 



This valve is a typical V-ported positioner equipped 
diaphragm control valve, which has as input the control signal Pdw 
and according to this signal the valve controls the amount of feedwater 
flow, Gw. The time delay shown represents the sticking of the valve 
stem . 

3 . Water Flow Transmitter 

This transmitter measures the differential pressure related 
to the flow of feedwater, Gw, across an orifice and develops a pneu- 
matic loading pressure Pgw, which is proportional to the feedwater flow. 

4 . Water Flow Feedback Signal Filter 

This is a needle valve-volume tank that serves the function 
of an R-C filter in the feedwater flow feedback loop, and is designed 
to attenuate the control system response to high-frequency components 
in the measured water flow signal, Pgs. The fact that this filter 
appears in the feedback portion of the loop causes a loop response 
similar to that obtained with proportional plus derivative response in 
the forward path. This in turn means that adjustment of the filter 
natural frequency will produce significant effects on the dynamic 
characteristics of the closed loop. 

5 . Steam Flow - Water Flow Differential Relay 

This relay subtracts the pneumatic signal, Psw', output 
of the filter, from the signal Pgs, output of the steam flow transmitter. 
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and develops a pneumatic signal, Psw’ , which is linearly proportional 



to this difference; the output of the relay is sent to the water level 
controller . 

B. ANALYSIS AND SIMULATION 

Figure 17 is the step response of the water loop with the param- 
eters and transfer functions as given in the block diagram. Figure 16. 

As before a parameter plane study was made to attempt the improvement 
of the loop step response. 

Figure 18 is the parameter plane plot for the loop made by 
defining the gains of the controller as: 

LAMBDA + RHO/S 

where , as given in Figure 16 , Lambda = 1.26, Rho = 0.10, and 
operating point marked "B" for which zeta = 0.21, Wn = 0.828 and 
sigma = 0.033. The parameter plane for this loop has certain resem- 
blance with the one for the oil loop, with the major difference that the 
complex roots are much closer to the origin in the present case making 
the loop more oscillatory. Then to damp the oscillation, zeta should 
increase or Lambda decrease which will also increase Wn; by increasing 
Rho, sigma is increased. The overall result of these adjustments is 
the transient response shown in Figure 19 made for Lambda = 0.5, 

Rho = 0.2, and it is seen that the response is faster and less 
oscilla tory . 
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WATER FLOW IN LB/HR X 1000 




FIGURE 17. UNIT STEP RESPONSE 
WATER LEVEL LOOP. LAMBDA = 1.26, RHO = 0.10 
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FIGURE 18. PARAMETER PLANE 
WATER LEVEL CONTROL LOOP 
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WATER FLOW LB/HR X 1000 
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Figure 19 

UNIT STEP RESPONSE WATER LEVEL CONTROL LOOP 
LAMBDA = 0.50, RHO = 0.20 
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V. SIMULATION OF COMPLETE SYSTEM 



A. DESCRIPTION 

The complete boiler control system consists, in addition to the 
already described components, of the following ones: 

1 . Steam Flow Transmitter 

This transmitter is identical to that of the water flow. The 
input is the load, Gs, steam flow; its output is proportional to the 
steam flow and is sent to the water level control loop and the steam 
pressure controller. 

2 . Steam Pressure Controller 

This instrument is a summing device that performs the 
operation KPep + Pgs, and the output is the control signal Pm, that 
is the input to the air loop, which is sent also to the signal selector 
relay. 

3 . Signal Selector Relay 

This relay selects the smaller of the two signals: Pm, 
from the steam pressure controller, and Pbg from the air loop, and 
sends it as demand index, to the fuel oil control loop. The reason for 
this is that the quantity of fuel oil supplied to the boiler furnace, in 
this way, can never be excessive for the amount of air flow measured 
by the air flow transmitter. In this manner, smoking is presented 
under conditions of rapidly increasing boiler load. 
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4 . Main Propulsion Steam Generator 

The steam generator was modeled by the NBTL [1] by 
simulating triangular pulses for the steam flow, the fuel oil flow and 
the feedwater flow, and measuring the response of steam drum pressure 
and water level to these forcing triangular pulses. As a result the 
steam generator was represented by a total of seven parallel transfer 
functions, three representing the dynamic behaviour of the steam 
pressure due to disturbances in steam flow, fuel oil flow and feedwater 
flow respectively, and four representing the dynamic behaviour of the 
water level to the same disturbances. As explained in Ref. 1, the 
response of the water level due to steam flow disturbances was sepa- 
rated in two parallel transfer functions, one attributable to the mass 
balance integration, since the boiler water level integrates at a rate 
proportional to the steam flow - water flow difference, the other 
attributable to the boiler water level "shrink and swell" . 

5 . Water Level Transmitter 

This a bellows type inverse acting differential pressure 
transmitter. It is designed to measure and transmit differential head 
pressure between water level in the boiler drum and a reference column 
of condensate connected into the steam space above the water level in 
the steam drum; the steady state calibration curve changes one psi for 
each one inch of boiler level change. The output, PI, is sent to the 
water level controller. 
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6 . Superheater Pressure Transmitter 



Measures the superheater outlet pressure and produces 
an output, Pp, between 3 and 2 7 psi which is compared against the 
reference signal Prp, set at 15 psi, and the difference is the input to 
the steam pressure controller. 

7 . Superheater 

The pressure drop across the superheater was represented 
as a linear function, based on the derivatives of the pressure drop 
curves evaluated at the two load conditions. In reality these pressure 
differentials vary with the square of the steam flow, and this fact is 
taken into account for the simulation of part VI of this thesis. 

B. SIMULATION 

The system was modeled by NBTL, as illustrated in Figure 1, 
at two load conditions, the cruising condition (Gs = 56000 Ib/hr) and 
the 90% of full power (Gs = 152000 Ibs/hr) . No data was found con- 
cerning the behaviour of the actual boiler to small perturbations of the 
load around the two conditions given, to be compared with the results 
of this simulation; therefore, a forcing function was chosen similar 
to the one used in Ref. 5 so that a comparison could be made, at least, 
of two different simulation methods for the same system. 

1 . Increasing Steam Load Condition 

For the condition of increasing load, a positive ramp was 
used for Gs. When the steam load on the boiler increases, the initial 
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effect is a decrease in superheater outlet pressure. As the pressure 
decreases, the output signal of the steam pressure transmitter decreases, 
causing a corresponding increase in the output signal from the steam 
pressure controller. 

The increasing signal from the steam pressure controller 
represents an increased signal of air demand. Since the air demand 
signal increases, the signal of air flow is less than the signal of air 
demand, and the air flow controller transmits an increased correcting 
signal. The increased signal pressure, in turn, increases blower 
speed and opens the forced draft dampers to increase air flow across 
the boiler air registers. 

On increasing the boiler load, the signal of air flow is 
less than the signal of boiler load; therefore, the signal of air flow 
is transmitted by the minimum signal selector relay as the minimum of 
the two signal pressures. 

As air flow across the air registers increases, the in- 
creased flow is sensed by the air flow transmitter which sends an 
increasing signal through the control circuit to increase air loading to 
the fuel oil pressure control valve. The increased loading on the 
diaphragm of the fuel oil control valve closes off the valve to increase 
the oil pressure in the fuel oil return line which increases the flow of 
oil to the burners to the amount required for optimum combustion. 

When the firing rate of the boiler has been increased to 
the level required to restore superheater outlet pressure to set point. 
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the combustion control system remains at the steady-state condition 



until again disturbed by variations in boiler load. 

a. Cruising Condition Response 

The load perturbation used for this condition and 
shown in Figure 20A was as follows: the steam flow was set initially 
at 56000 Ib/hr and held for 60 seconds, so that all variables of the 
system reach steady-state condition; then it was ramped to 59000 Ib/hr 
in five seconds and held at that value for fifty seconds. The responses 
were recorded as follows: 

(1) Water Level Response . Figure 20B, where it 
is seen that a peak of 0.19 inches is reached in 19 seconds and then 
the level approaches zero again. 

(2) Superheater Outlet Pressure Response . 

Figure 20C. It has a minimum of -4.32 psi reached in 10 seconds 
after the perturbation was applied and then it approaches the steady 
state which is seen to be about one psi above 1200 psi. 

(3) Air Flow Response . Figure 2 0D. It goes from 
19.1% to 20.6% following the ramp with an overshoot that has a peak 
of 21% at 7 seconds, and then approaches the steady state. 

(4) Fuel Oil Flow Response . Figure 2 0E. This 
response is similar in shape to the air flow response, has an overshoot 
that peaks at 17 seconds with a value of 4400 Ib/hr and then settles to 
a value of 4230 Ib/hr in a much faster rate than the air flow does. 
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STEAM PRESSURE WATER LEVEL STEAM FLOW 




10 20 30 40 50 

Figure 20A. Increasing Steam Load (Cruising) 

Given in Ib/hr x 1000.0, vs. time in seconds 



CM 




Figure 20B. Water Level Response 
Given in inches of variation from set level, vs. time in seconds 




Figure 20C. Steam Pressure Response 
Given in psi variation from 1200psi, vs. time in seconds 
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FEEDWATER FLOW FUEL OIL FLOW AIR FLOW 




Figure 20D. Air Flow Response 
Given in percentage of full power, vs. time in seconds 




Figure 20E. Fuel Oil Flow Response, vs. time in seconds 
Given in Ib/hr variation from steady state value (4040.) 




Figure 20F. Feedwater Flow Response, vs. time in seconds 
Given in Ib/hr x 1000.0 variation from 56000.0 Ib/hr 
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(5) Feedwater Flow Response . Figure 20F. The 



response of the feedwater flow is seen to be very slow, and oscillatory 
at the beginning of the transient. 

All of the responses specified above are in agreement 
with the type of response expected from the discussion given at the 
beginning of this part. If these responses are compared with the 
corresponding ones given in Ref. 5, it is seen that they are very 
similar; the shapes of the responses are the same. The major difference 
is noted in the feedwater response, which for the case of this simulation 
is much less oscillatory, due to the changes made in the transfer 
functions . 

b. 90% Full Power Condition Responses 

The load perturbation used for this condition is shown 
in Figure 21A and was simulated as follows; after steady state was 
reached for 152000 Ibs/hr of steam flow it was ramped to 160000 Ibs/hr 
in five seconds and held for fifty seconds. The corresponding responses 
were: 

(1) Water Level Response . Figure 2 IB. It shows 

a peak of 0.38 inches at 16 seconds after the perturbation was applied, 
and then the steady state is approached. 

(2) Superheater Outlet Pressure Response . 

Figure 2 1C. It has a minimum of -3.8 psi at 8 seconds, and then goes 
to the steady state in about 9 seconds more, with a small error of less 
than one psl above 1200 psi. 
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STEAM PRESSURE WATER LEVEL STEAM FLOW 




Figure 21A. Increasing steam load (90% FP) 
Given in Ib/hr x 1000.0, vs. time in seconds 




Figure 2 IB. Water Level Response, vs. time in seconds 
Given in inches of variation from set level 




Figure 2 1C. Steam Pressure Response, vs. time in seconds 
Given in psi variation from steady state 
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FEEDWATER FLOW FUEL OIL FLOW AIR FLOW 



oi 

U3 




10 20 30 40 50 



Figure 2 ID. Air Flow Response, vs. time in seconds 
Given in percentage of full power 




Figure 2 IE. Fuel Oil Response, vs. time in seconds 
Given in Ib/hr from steady state (124 80 Ib/hr) 




Figure 2 IF. Feedwater Flow Response, vs. time in seconds 
Given in Ib/hr x 1000 variation from 152000 Ib/hr) 
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(3) Air Flow Response « Figure 2 ID. It goes from 



59% to 62% in 2 5 seconds, with an overshoot of 64.2% at 12 seconds. 

(4) Fuel Oil Flow Response . Figure 2 IE. Goes 
from 12480 Ib/hr of fuel oil flow to 13100 Ib/hr in about 23 seconds, 
overshooting to 13440 at 12 seconds. 

(5) Feedwater Flow Response . Figure 2 IF. Very 
slow to reach steady state and also shows some oscillation at the 
beginning of the transient. 

Comparison of these responses with the corresponding 
ones of Ref. 5 shows good agreement, with the following noticeable 
differences: variation in steam pressure is less for this simulation; 
the feedwater flow is also seen to be much less oscillatory for this 
simulation . 

As was said before these differences are attributable 
to the corrections made to some of the transfer functions. 

2 . Decreasing Steam Load Condition 

For this condition a negative ramp was used to simulate 
the load, Gs . When this reduction cf steam load occurs an instan- 
taneous increase in superheater outlet pressure takes place. With 
increasing pressure there is an increasing output signal from the steam 
pressure transmitter which causes a proportional decreasing signal 
from the steam pressure controller. This decreasing signal indicates 
a reduced signal of air demand at the air flow controller. Therefore, 
the signal of air flow will be greater than the signal of air demand 
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which causes a correcting decreased signal to be transmitted by the air 



flow controller. This decreased signal pressure decreases blower speed 
and closes off the forced draft dampers to reduce combustion air flow 
to the boilers . 

With decreasing boiler load, the boiler load signal to the 
minimum signal selector relay is instantaneously less than the signal 
of air flow. 

When the boiler load signal is less than the signal of air 
flow, the minimum signal selector relay immediately transmits the lower 
boiler load signal to reduce oil flow to the burners. 

As the firing rate is adjusted to restore superheater outlet 
pressure to set point, the automatic combustion control system returns 
to its steady-state condition until variations in boiler load again cause 
a change in superheater outlet pressure. 

a. Cruising Condition Response 

The load perturbation used for this condition is shown 
in Figure 22k, and was as follows: after steady state was reached by 
the system at 56000 Ib/hr of steam load a negative ramp was applied 
for five seconds to bring Gs to 53000 Ib/hr. The responses were as 
follows: 

(1) Water Level Response . Figure 22B. Has a 
maximum decrease of 0.185 inches at 11 seconds after the perturbation 
and then goes to the set level again. 
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Figure 22A. Decreasing Steam Load (cruising) 
Given in Ib/hr x 1000, vs. time in seconds 




Figure 22B. Water Level Response, vs. time in seconds 
Given in inches of variation from set level 




Figure 22C. Steam Pressure Response, vs. time in seconds 
Given in psi variation from steady state 
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FEEDWATER FLOW FUEL OIL FLOW AIR FLOW 




Figure 2 2D. Air Flow Response, vs. time in seconds 
Given in percentage of full power 
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Figure 22E. Fuel Oil Response, vs. time in seconds 
Given in Ib/hr variation from steady state (4040 Ib/hr) 




Figure 22F. Feedwater Flow Response, vs. time in seconds 
Given in Ib/hr x 1000 variation from 56000 Ib/hr 
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(2) Steam Pressure Response . Figure 22C. 



Initially the pressure increases to about 1 psi above the steady state 
in five seconds, then decreases to -1 .25 psi below the steady state 
and keeps this value as new steady state. 

(3) Air Flow Response . Figure 22D. Varies from 
19% to 18.2% in fifteen seconds. 

(4) Fuel Oil Flow Response . Figure 22E. Decreases 
from steady state with a small overshoot at seven seconds after the 
perturbation and then settles to a value of 3 860 Ib/hr. 

(5) Feedwater Flow Response . Figure 22F. Presents 
some oscillation during the first fifteen seconds of the response and then 
decreases steadily to its new steady state value of 53000 Ib/hr. 

b. 90% Full Power Condition Responses 

The load perturbation used, Figure 23A, was, as 
before, letting the system reach steady state at 152000 Ib/hr and then 
a negative ramp was applied for five seconds to bring the load, Gs, to 
144000 Ib/hr. The responses were as follows; 

(1) Water Level Response . Figure 23B. Shows 
a minimum of -0.3 9 inches at 11 seconds and then goes to the set 
level again. 

(2) Steam Pressure Response . Figure 23C. 

Initially the pressure increases to about 4 psi above steady-state 
value at 6 seconds, then decreases to its steady-state value in about 
six seconds more and settles at a value of about one psi below 1200 psi. 
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Figure 23A. Decreasing Steam Load (90% FP) 
Given in Ib/hr x 1000.0, vs. time in seconds 
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gure 23B. Water Level Response, vs. time in seconds 
Given in inches of variation from set level 




Figure 23C. Steam Pressure Response, vs. time in seconds 
Given in psi variation from steady state 
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Figure 23D. Air Flow Response, vs . time in seconds 
Given in percentage of full power 
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Figure 23E. Fuel Oil Flow Response, vs. time in seconds 
Given in Ib/hr variation from steady state value (12480 Ib/hr) 




Figure 23F. Feedwater Flow Response, vs. time in seconds 
Given in Ib/hr x 1000 variation from 152000 Ib/hr 
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(3) Air Flow Response . Figure 23D. Goes in eight 



seconds from 59.0% to 55.7%, its new steady state. 

(4) Fuel Oil Flow Response . Figure 23E. 

Decreases from steady-state value (12480 Ib/hr) to its new steady 
state (11790 Ib/hr) in about 10 seconds. 

(5) Feedwater Flow Response . Figure 23F. Again 
shows oscillation during the beginning of the response and then decreases 
steadily to its new value of 144000 Ibs/hr. 

Comparison of the decreasing responses to the 
corresponding ones of Ref. 5 shows good agreement, with the same 
exception as before of the feedwater flow response. For the steam 
pressure response to a decreasing ramp, 90% FP condition, in Figure 
7-12 (a) of Ref. 5 the steam pressure at the beginning of the response 
is shown to decrease, which is not consistent with the physical 
behaviour of steam pressure under decreasing load condition. 

3 . New Parameters Simulation 

The simulation of the system using small perturbations 
was repeated using as gain settings for the different controllers, the 
values found in parts II, III, and IV of this thesis. The general shape 
and values of the responses were very close to those obtained using 
the original gain settings. The most noticeable change for all of the 
responses is that the feedwater response was less oscillatory and 
somewhat faster. 
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Increasing Steam Load 



a . 

(1) Cruising Condition . Shown in Figures 24A 
to 24F. Responses are very close to those previously obtained; the 
feedwater response is less oscillatory. 

(2) 90% Full Power Condition . Figures 25A to 
25F. Same behaviour as noted for cruising condition; here the feed- 
water response is noticeable faster than before. 

b. Decreasing Steam Load 

(1) Cruising Condition . Figures 26A to 2 6F . Same 
observations as noted for the increasing load condition. 

(2) 90% Full Power Condition . Figures 2 7Ato 
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Figure 24A. Increasing Steam Load (cruising) 
Given in Ib/hr x 1000.0, vs. time in seconds 




Figure 24B. Water Level Response, vs. time in seconds 
Given in inches of variation from set level 




Figure 24C. Steam Pressure Response, vs. time in seconds 
Given in psi variation from steady state 
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Figure 24D. Air Flow Response, vs. time in seconds 
Given in percentage of full power 




Figure 24E. Fuel Oil Flow Response, vs. time in seconds 
Given in Ib/hr variation from steady state 




Figure 24F. Feedwater Flow Response, vs. time in seconds 
Given in Ib/hr x 1000 variation from 5 6000 Ib/hr 
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Figure 25A. Increasing Steam Load (90% FP) 
Given in Ib/hr x 1000.0, vs. time in seconds 
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Figure 25B. Water Level Response, vs. time in seconds 
Given in inches of variation from set level 
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Figure 2 5C. Steam Pressure Response, vs. time in seconds 
Given in psi variation from steady state 
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uel Oil Flow Response, vs, time in seconds 



Given in Ib/hr x 1000 variation from steady state 




Figure 25F. Feedwater Flow Response, vs. time in seconds 
Given in Ib/hr x 1000 variation from 152000 Ib/hr 
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Figure 26A. Decreasing Steam Load (cruising) 
Given in Ib/hr x 1000.0, vs. time in seconds 




Figure 26B. Water Level Response, vs. time in seconds 
Given in inches of variation from set level 




Figure 2 6C. Steam Pressure Response, vs. time in seconds 
Given in psi variation from steady state 
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Figure 26E. Fuel Oil Flow Response, vs. time in seconds 
Given in Ib/hr variation from steady state 




Figure 26F. Feedwater Flow Response, vs. time in seconds 
Given in Ib/hr x 1000 variation from 56000 Ib/hr 
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STEAM PRESSURE WATER LEVEL STEAM FLOW 




Figure 27A. Decreasing Steam Load (90% FP) 
Given in Ib/hr x 1000.0, vs. time in seconds 




Figure 27B. Water Level Response, vs. time in seconds 
Given in inches of variation from set level 




10 20 30 40 50 



Figure 2 7C. Steam Pressure Response, vs. time in seconds 
Given in spi variation from steady state 
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Figure 27D. Air Flow Response , vs. time in seconds 
Given in percentage of full power 




Figure 27E. Fuel Oil Flow Response, vs. time in seconds 
Given in Ib/hr x 1000 variation from steady state 




Figure 27F. Feedwater Flow Response, vs. time in seconds 
Given in Ib/hr x 1000 variation from 152000 Ib/hr 
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VI. LINEAR INTERPOLATION SIMULATION 



A. DISCUSSION 

As already has been discussed the nonlinearities of the model are: 
-The gain of the blower actuators 
-The dynamics of the main forced draft blowers 
-The seven transfer functions representing the steam 
generator 

-The pressure drop across the superheater 
The pressure drop across the superheater is represented as a 
linear function in the model, but this pressure drop varies with the 
square of the steam flow, as explained before. The following expres- 
sion was used in the simulation to represent this pressure drop; 

Pp = (Gs)^/(1 .86x10®) 

where Pp is the superheater outlet pressure in psi, and Gs is the steam 
flow in Ib/hr. The constant of proportionality was chosen as the 

g 

inverse of 1.86 x 10 , so that the two points given in the model for 

cruising and 90% FP are contained by the square law curve; that is: 

-4 2 

for cruising: (56000.0) x (3x10 ) = (56000.) x K 

for 90% FP: (152000. )x(8 .ZxlO”"^) = (152000.)^ x K 

g 

Solving for K gives: for cruising K = 1.87x10 and for 90% FP K = 
1.85x10®. 

All the other gains and time constants of the nonlinear functions 
were simulated as follows: 



82 



For loads of 90% FP and above, the values given for the 90% FP 



condition were used. 

For loads between 90% FP and Cruising conditions, a linear 
interpolation was made between the values given for these two con- 
ditions by the NBTL report. 

For loads below Cruising (56000. Ib/hr) an extrapolation was 
made continuing the same straight line used for the interpolation, 
down to a value of Gs equal to 28000. Ib/hr. 

The extrapolation was not continued for lower values of Gs , or 
above 90% FP, because some values get very close to the imaginary 
axis and the results obtained were very oscillatory. 

For loads below 28000. Ib/hr the values of gains and time con- 
stants were held constants at the values calculated for 28000.0 Ib/hr. 

B . SIMULATION 

The simulation was done using a forcing function that repre- 
sented a linear variation of load from 10% FP to 90% FP and vice versa, 
which was the same load variation used in Ref. 4 to drive a test DLG-9 
Class boiler. 

For comparison, the responses to the same driving force, but 
using the cruising condition transfer functions only, and the 90% FP 
transfer functions only, are also shown, along with the interpolation 
and the actual responses. 
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1 . 



Increasing Load Condition 



As before, the system was allowed to reach steady-state 
condition for all variables by applying a load of 10% FP, that is, 

16880. Ib/hr of steam flow, for a period of 60 seconds. Then a linear 
increasing ramp was applied up to 90% FP, in 23 seconds. The responses 
to this load variation are shown in Figures 28 to 31, where the first 
graph of each figure is the response of the test boiler to the same load 
disturbance, as reported in Ref. 4. The second graph is the response 
using for the model the cruising condition transfer functions only; the 
third graph is the response of the model using the 90% FP transfer 
functions only and the fourth graph is the response of the model using 
the interpolation as explained above. 

Water Level response, (Figure 28). The model responses 
are similar in shape to the actual response; they peak at about 28 
seconds, and then tend to the normal level. However, all of the model 
responses are larger in magnitude than the actual response. 

Steam Pressure response, (Figure 29). The actual 
response peaks at about 23 seconds, with a value of 1115 psi while 
the model responses peak at about 16 seconds with the following 
magnitudes: Cruising, 1125 psi; 90% FP, 1144 psi; and interpolation, 
1082 psi. The actual response settles to 1200 psi, while the model 
responses settle to values higher than 1200 psi, as seen in Figure 29. 

Air Flow response, (Figure 30). Model responses are 

similar in shape to the actual, but cruising and 90% FP responses are 
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Figure 28. Water Level Response to increasing load 
Inches of variation from set level N, vs. time in seconds 
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Figure 29. Steam pressure response to increasing load 
Given in psi, vs. time in seconds 
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Figure 30. Air Flow response to increasing load 
Given in percentage of FP, vs. time in seconds 
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Figure 31. Oil Flow response to increasing load 
Given in percentage of FP, vs. time in seconds 
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Figure 32 , Water Flow response to increasing load 
Given in percentage of FP, vs. time in seconds 
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smaller in magnitude of peak and settling values. The interpolation 
response values are closer to the actual response values than the 
cruising and 90% FP values are. 

Fuel Oil Flow response, (Figure 31). The observations 
made for the air flow response are valid also for the fuel oil flow 
response. Interpolation responses are closer to actual than the other 
two simulations . 

Water Flow response , (Figure 32) . The actual response 
is observed to begin at about 13 seconds while the model simulations 
begin responses at about 4 seconds; also the actual response settles 
to 60% at about 50 seconds, while the model responses are slower to 
reach the steady state . 

2 . Decreasing Load Condition 

The system was driven with a steady load of 152000. Ib/hr 
of steam load for 60 seconds until steady-state values were reached 
for variables; then a steady decreasing ramp was applied down to 10% 
of full power during 23 seconds. Figures 33 to 37 show the responses, 
in the same form as for the increasing load condition. 

Water level response, (Figure 33). Actual response of 
water level was a decreasing of level of 3.9 inches in 25 seconds; 
the model responses show a decreasing of level of more than 5.5 inches 
in 2 7 seconds for all the three simulations. 

Steam pressure response, (Figure 34). Actual response 
shows that superheater pressure goes from steady state (1210 psi) to 
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Figure 33. Water level response to decreasing load 
Inches of variation from set level N, vs. time in seconds 
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Figure 34. Steam pressure response to decreasing load 
Given in psi, vs. time in seconds 
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Figure 35. Air Flow response to decreasing load 
Given in percentage of FP, vs. time in seconds 
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Figure 36. Oil flow response to decreasing load 
Given in percentage of FP, vs. time in seconds 
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Figure 37. Water flow response to decreasing load 
Given in percentage of FP, vs. time in seconds 
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1248 psi in 8 seconds and reaches a new steady state of 1185 psi at 
40 seconds. Model simulations were as follows; Cruising, from 1235 
psi to 1248 psi in 7 seconds, and settles to 1195 psi at 28 seconds; 

90% FP, from 1220 psi to 1230 psi in 7 seconds, and reaches steady 
state of 1200 psi at 35 seconds with small oscillations around this 
value; interpolation response was from 1215 psi to 1310 psi in 10 
seconds, reaching new steady state of 1200 psi at 40 seconds. 

Air flow response, (Figure 35). Responses of the model 
were the same as actual in shape but magnitudes of the model were 
smaller than the actual. 

Oil flow response (Figure 36). Again the magnitudes of 
model responses were smaller than magnitudes of actual responses. 

Water flow response, (Figure 37). Same observation as 
for the increasing load condition. The water flow responses were 
faster to initiate their response, but slower to reach the steady state 
than the actual responses. 

C. NEW PARAMETERS SIMULATION 

Comparison of the responses presented in the preceding part B, 
shows that the model with 90% FP transfer functions is a better model 
to represent the steam generator and its control system, for large 
variations of load, than the cruising condition or the linear interpolation 
model. Therefore, the 90% FP transfer functions were used to perform 
a simulation using the new parameters and variations of load from 10% 
to 90% of full power, and vice versa. 
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Of course, there are no actual boiler responses for the new param- 



eters to compare with; such is the case of the designer who is presented 
with a model of the plant and he is to find the controller gains, to make 
the plant respond according to certain specifications. 

In the discussion to follow, comparison will be made with the 
90% FP responses of the model using the original parameters. 

1 . Increasing Load Condition 

The forcing function is the same as the one used in part 
B, that is, the system is allowed to reach steady-state values at 10% 
of full power, and then it is ramped linearly, and in 23 seconds, to 
90% of full power, as shown in Figure 38A. 

The responses, shown in Figures 38B to 38F, are not 
significantly different from those for the original parameters, they are 
somewhat faster in the sense that the peaks and settling times are 
about 3 seconds earlier than before. The steam pressure response 
peaks about 6 seconds earlier than previously. 

2 . Decreasing Steam Load 

The forcing function is shown in Figure 39A and the 
responses are shown in Figures 39B to 39F; they are similar in shape, 
magnitude and times of peaking to the responses using the old param- 
eters, however, the settling time is increased since some oscillation 
is noted at the latter part of the response. 
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Figure 38A. Increasing Steam Load 
Given in percentage of FP, vs. time in seconds 
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Figure 38B. Water Level Response 
Inches of variation from set level N, vs. time in seconds 




Figure 38C. Steam Pressure Response 
Given in psi , vs. time in seconds 
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Figure 38D. Air Flow Response 
Given in percentage of FP, vs. time in seconds 




Figure 38E. Fuel Oil Flow Response 
Given in percentage of FP, vs. time in seconds 




Figure 38F. Water Flow Response 
Given in percentage of FP, vs. time in seconds 
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Figure 39A. Decreasing Steam Load 
Given in percentage of FP, vs. time in seconds 




Figure 39B. Water Level Response 
Inches of variation from set level N, vs. time in seconds 




Figure 39C. Steam Pressure Response 
Given in psi, vs. time in seconds 
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Figure 39D. Air Flow Response 



Given in percentage of FP, vs. time in seconds 




Figure 39E. Fuel Oil Flow Response 
Given in percentage of FP, vs. time in seconds 




Figure 39F . Water Flow Response 
Given in percentage of FP, vs. time in seconds 
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VII. CONCLUSIONS 



A. DISCUSSION OF RESULTS 

Parts II, III, and IV of this thesis demonstrate the use of 
parameter plane techniques, as applied to complex control loops, 
where all the parameters of the loop had been selected and there 
remains only the problem of selecting the parameter settings (gain 
values) of the controllers to provide acceptable, (according to certain 
criterion), loop performance. 

The use of the parameter plane will reduce by a significant 
amount the time expended in field trials to find these gain settings. 

It is to be noted that the optimization of several subloops 
according to certain specifications, as was done in parts II, III, and 
IV does not imply that when these subloops are put together the entire 
system will perform as expected from the specifications used in the 
optimization of the subloops, because of the many interactions between 
them . 

Part V shows that the behaviour of the system to small dis- 
turbances was as expected from theory; however, the magnitude of 
the variations of the responses could not be compared with responses 
of actual boiler to small perturbations, because of lack of relevant 
data . 

Part VI assesses the validity of the model, since the responses 
of this mathematical model arc close to the responses of the actual 
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boiler. However, some discrepancies in magnitude were noted, because 



the transfer functions were found at two fixed operating conditions, and 
hence the model is valid only for small perturbations about these 
operating conditions. 

The responses of the model shown on part VI indicate that the 
90% Full Power condition transfer functions provide a better model for 
large excursions of load, than the cruising condition transfer functions. 
It also shows that the linear interpolation did not bring the magnitudes 
of the responses closer to actual values than those for the 90% FP 
transfer functions, and in the case of steam pressure the variations 
were larger than when using cruising transfer functions. 

The use of the DSL/3 60 digital simulation language was very 
useful, and easy to use, to obtain the time responses of the model. 

B. RECOMMENDATION 

In order to have a model that closely represents the behaviour 
of the actual boiler for large changes in steam flow, further work is 
recommended in the interpolation and extrapolation of the nonlinear 
transfer functions, since the linear approximation done in part VI of 
this thesis did not improve the 90% FP transfer functions responses. 
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APPENDIX A 



PERSONAL LETTER 




DEPARTMENT OF THE NAVY 



NAVAL SHIP ENGINEERING CENTER. PHILADELPHIA DIVISION 
PHILADELPHIA, PA. 19112 



IN REPLY REFER TO 



6770:JWB:jm 
Ser 30 
3 Feb 1970 



( 



LCDR Humberto Medina, Venezuelan Navy 
U. S. Naval Postgraduate School 
SMC No. 1486 
Monterey , Calif. 93940 

Dear LCDR Medina: 

Thank you for your interesting letter concerning your simulation 
studies of the DLG-9 Class steam generator. It appears that you have 
indeed identified one of the trouble spots In our "model" of this 
system. This same difficulty was not recognized by us until 1966 
when we operated a control system on the boiler; the system had been 
designed with the aid of simulation based upon the "model" reported in 
our Project No. B-502-III, 

Let us first clear up your confusion concerning our block diagram. 
Figure 3 of the report. It Is true that the constant term of the 
characteristic of the closed loop fuel system is negative; this is com- 
pensated by the fact that the fuel control valve, Gy, is reverse._ac.ting, 
thus placing a negative sign In front of the dynamic terms of both 
numerator and denominator of the closed loop transfer function. The 
negative signs of course cancel out, leaving us with a potentially 
stable system. You have correctly concluded that the system is not 
stable for the controller settings indicated. This results from the 
fact that the elements of the control loop were individually identified 
by open- loop testing; as a result, the effects of impedance mismatching 
on the system dynamics were not revealed. It turns out that the loading 
effect of the fuel control valve diaphragm operator on the fuel controller 
effectively adds a tv/o-second first order time constant to the system. 

Thus the control valve transfer function should be properly specified 
as Gv = -795/(1 + 2. OS). This problem is discussed in some detail In 
our Report No. B-584-II, which may be obtained from the Defense 
Document Center using AD No. 805-700L. I believe that you will find 
that the addition of this first order lag to the system will stabilize 
your simulation. You should not have to alter the dynamics of either 
the controller or the transmitter to achieve stability. 
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In addition to the correction to the high frequency components of 
the water level controller, the only other correction to the system is 
that of the gain of the feedwater control valve. It should be changed 
from 3.12 x 10** to 0.94 x 10**. The new figure represents the average 
gain of the valve over the full load range, whereas the original figure 
represents a perturbation result at one load level. 

The control devices in the system were manufactured by the Bailey 
Meter Company; all operate on a signal range of 3 to 27 psig. The 
controller settings were obtained by tria l-and-error adjustment in the 
field based on transient response specifications of +5 in. of water 
level and -120 psi superheater outlet pressure as load varies from 8 
to 15%, and -5 in., +20 psi as load varies from 75 to 8^, both in 
25 seconds and linearly. Remember that pressure drop through the super- 
heater varies with the square of steam flow and amounts to 75 psi 
at 83.3^ 9 f rating. 

I hope that the foregoing information will be of value to you in 
your v/ork; please feel free to correspond if additional assistance is 
requ i red. 



Very truly yours. 




Head, Machinery Automation 
Systems Department 



105 



APPENDIX B 



PERSONAL LETTER 



PHILADELPHIA NAVAL SHIPYARD 

PHILADELPHIA. PA. 19112 in reply refer to 

Code 2725 



30 August 1965 



Professor I.lilton L. V/ilcox 
U. S. Naval Postgraduate School 
Itonterey, Calif. 93940 

Dear Professor Wilcox: 

Thank you for your letter of 10 August 1965 expressing your interest in 
a continuation of the simulation program undertaken by Lieutenants Creager, 
Fenickj and O'Brien. I apologize for not having replied sooner, but I was 
on my vacation v/hen your letter arrived. 

I v/ould be very pleased to provide whatever assistance and information 
required that is available here at the Laboratory. In reply to your 
specific questions, the follov/ing information is offered: 

a. There is no additional test information available on the 
feedwater system. 

b. The entire boiler and control system has been satisfactorily 
simulated on the analog computer here at NBTL. 

c. There are errors in the transfer functions shown in blocks 29 
and 84 of Figure 6-3 of the simulation thesis. These errors resulted 
from erroneous information originally provided by xis, but I had been 
under the impression that it'v/as subsequently corrected. The correct 
transfer functions should be 12 5/ (S^ + 6.75S + 125) for block 84, and 
19.2/(s2 + 2.62S + 19.2) for block 29. These errors undoubtedly account 
for the instability of the simulated digital system. It should be noted 
that the analog simulation v/as stable, and, in general, duplicated the 
transient response of the actual boiler. 

V/e are presently preparing a detailed description of the procedure whereby 
the analog program was developed. A copy of the computer program of the 
feedv/ater control system, together vn.th this description, will be forwarded 
to you as soon as completed. 
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I have been in touch with };tr. Lackowski concerning our mutual problems 
and interests-; subject to authorization by the Bureau of Ships, some 
sort of joint effort betv/een our tv/o Laboratories may result from this 
contact. 

As your work progresses, should you find need for additional information, 
please advise LCDR Ediin and LT Lamb to feel free to contact us directly. 

I suggest that any correspondence be addressed officially to the Commander, 
Philadelphia Naval Shipyard (Naval Boiler and Turbine Laboratory), and 
marked to ray attention. In this manner we can maintain such correspon- 
dence as a matter of record. 



Very truly yours, 




L. 



c JA'ffiS W. BANMi, JR. 
Head, Controls Branch 
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